Macroscopic Quantum Tunneling in Small Antiferromagnetic Particles: Effects of a 

Strong Magnetic Field 
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We consider an effect of a strong magnetic fieid on the ground state and macroscopic coherent 
tunneling in small antiferromagnetic particles with uniaxial and biaxial single-ion anisotropy. We 
find several tunneling regimes that depend on the direction of the magnetic field with respect to 
the anisotropy axes. For the case of a purely uniaxial symmetry and the field directed along the 
easy axis, an exact instanton solution with two different scales in imaginary time is constructed. 
For a rhombic anisotropy the effect of the field strongly depends on its orientation: with the field 
increasing, the tunneling rate increases or decreases for the field parallel to the easy or medium axis, 
respectively. The analytic results are complemented by numerical simulations. 
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During the last decade macroscopic coherent quantum 
tunneling between two equivalent classical ground states 
in macroscopic or, to be more precise, mesoscopic mag- 
netic systems became an object of intense experimental 
and theoretical investigations, see for a review Refsii^. 
In the physics of magnetism such systems are, for in- 
stance, small magnetic particles, magnetic clusters, and 
high-spin molecules. The interest in tunneling effects 
is associated with the possibility of using these systems 
as potential elements for quantum computers. Initially, 
the calculations of tunneling effects were carried out for 
ferromagnetsi&i However it happened that the tunneling 
effects were experimentally observed by the resonant ab- 
sorption of electromagnetic waves for antiferromagnetic 
ferritin particles^ The effect of the magnetic field on 
the tunneling probability in ferritin particles have been 
experimentally studied, see Ref£ According to theoreti- 
cal estimations the level spitting in antiferromagnets 
is stronger than in ferromagnets and the effects can be 
observed at higher temperatures. Antiferromagnets are 
more convenient for experimental investigations of tun- 
neling effects. 

The interest in tunneling in compensated antiferro- 
magnets has been renewed after the synthesis of high- 
spin molecules with antiferromagnetic coupling of spins, 
so-called ferric wheels, such as Fe6, Feio, Crg, see, for in- 
stance, RefsAiiii The interference effects caused by the 
gyroscopic terms in the a-model approximation were of 
main interest^ For compensated antiferromagnets such 
effects can arise due to an external magnetic fieldiS*^ or 
a certain kind of the Dzyaloshinskii-Moriya interaction^ 
Chiolero and Loss^ have also found a type of interfer- 
ence effect for such antiferromagnets placed in a strong 
magnetic field directed perpendicular to the easy axis of 
the antiferromagnetic particle. This effect caused by the 
imaginary part of the fluctuation determinant together 
with the direct contribution of the field to the imaginary 
part of the Euclidean action produces the periodic depen- 
dence of the level splitting on the field strength. Then Hu 
et al^ have investigated these effects for a more general 
case of antiferromagnets with rhombic anisotropy. They 



have considered the case of the field directed along the 
hard axis and found the behavior common to that in 
Ref^. They also suggested that such properties are uni- 
versal for a wide class of antiferromagnets subject to a 
strong magnetic field. 

In this paper we consider the effects of coherent quan- 
tum tunneling for a compensated antiferromagnetic par- 
ticle with rhombic anisotropy placed in the magnetic 
field. In the section [I] we describe the model used for 
the evaluation of tunnel splitting and analyze the static 
energy of the particle for three different orientations of 
the magnetic field along the crystalline axes. The sec- 
tion [n] is devoted to the calculation of the level splitting 
of the ground states in the instanton approximation for 
the orientations of the field along the easy, middle and 
hard axis separately and for relevant ranges of the field 
strength. The section [TTTl contains a qualitative compar- 
ison of analytical results with data of numerical calcula- 
tions for level splitting and a direct comparison for the 
Euclidean action. The final section [Tvl resumes obtained 
results and gives a short survey of experiments where the 
predicted effects could be observed. We use the standard 
semiclassical analysis based on the instanton technique 
applied to the nonlinear er-model, as well as the direct 
numerical diagonalization of the corresponding quantum 
spin Hamiltonian. We demonstrate that the behavior of 
the level splitting is highly sensitive to the orientation of 
the magnetic field. In fact, three different scenarios for 
the orientation of the field along the symmetry axes are 
found. 



I. NONLINEAR a-MODEL FOR 
ANTIFERROMAGNETS IN A MAGNETIC 
FIELD 

We start from the Hamiltonian describing a magnetic 
particle with an even number N of magnetic ions with 
the spin S and coupled with the nearest-neighbor antifer- 
romagnetic exchange interaction J. We also assume that 
a single-ion anisotropy with rhombic symmetry and the 
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magnetic field H are present. The macroscopic Hamil- 
tonian of the system can be written as 



<af}> 



+ B P Y,(S y a ) 2 -9VBj2 H - S *- (!) 



Here, S a is the spin at the ath site, g ss 2 is the Lande 
factor, and \ib is the Bohr magneton. The first term 
describes the isotropic exchange interaction J > 0, and 
the summation in this term is extended over the pairs 
of nearest neighbors. The second and third terms give 
the simplest form of a rhombic single-ion anisotropy; B u 
and Bp are the constants of the uniaxial anisotropy and 
the anisotropy in the basal plane xy, respectively. We 
assume that B u > and B p > 0. Thus, z is the easy axis 
and x, y are the medium and hard axes, respectively. 
The last term corresponds to the coupling spins with the 
magnetic field. 

We will analyze the system under the assumption that 
all spin pairs are equivalent and have the same coordina- 
tion number z. It is true for spin dimers (z = 1), spin 
wheels (z = 2), and can be a good approximation for 
mesoscopic three-dimensional crystalline particles with 
N 3> 1 , in which the surface variation of parameters can 
be neglected. Assume that the Zeeman energy g^igHS, 
the anisotropy energy B P S 2 , and B U S 2 are much smaller 
than the exchange energy JzS 2 . The classical approxi- 
mation gives that the spins in each sublattice are parallel. 
This assumption is a necessary condition for using the a- 
model as a model for describing both classic and quantum 
dynamics of antiferromagnets. 13 For the semiclassical dy- 
namics, especially for tunneling, it is convenient to use 
the Euclidean formulation, which is based on the intro- 
duction of imaginary time r = it. 

The Euclidean action for the cr-model that corresponds 
to the microscopic Hamiltonian can be written as 



/-t-oo 
C e [6{t), 0(t)] dr - 
-oo 

i 2 - 2i<yH ■ (I x i) 



N 



\4Jz 



W a (l) 



(2) 



where 7 = g/iB/h is the gyromagnetic ratio, N is the 
number of magnetic ions in both sublattices and the over- 
dot denotes the derivative with respect to the imaginary 
time r. 

In this approach the total spin Stot of the particle be- 
comes a slave variable, and it is determined through the 
unit Neel vector I and its derivative I with respect to r: 

Sua = ^ {7 [H - l(H ■l)]+i(lx /) } . (3) 

The function w a (l) is an effective energy of anisotropy 
per spin with renormalization caused by the external field 



w a (l) = B u S 2 (l 2 x + l 2 y ) + B p S 2 l 2 y 



AJz 



(H-lf. (4) 



In accordance with the general rules of the semiclassi- 
cal approximation formulated in the instanton language 
the amplitude of the tunnel transition from the state 
I = I to I = l + is proportional to exp(— Ae/K), where 
the value of Ae is calculated from appropriate equations 
of motion under the boundary conditions: i(r) — * I at 
r — > —00 and i(r) — > l + at r — > +00. Tunnel splitting 
in the so-called dilute instanton gas approximation is de- 
termined by the modulus of the sum of such amplitudes 
calculated along all instanton trajectories with a minimal 
value of Re Ae- 

It is convenient to introduce a polar parametrization 
for the unit vector I 



l z = cosO, l x = sin#cos(/>, l y — sin 9 sine 



(5) 



The Lagrangian from Eq. J2J in such a parametrization 
takes the form 



£ e [0(t),0(t)] =Nw a {9,4>) 

( ,0 y . 2 



4Jz 



+ </) sin 9 — 2ij < 9{H y cos 4> — H x sin < 



[H z sin 2 9 — (H x cos <ft + H y sin (f>) sin 9 cos 0] 1 



(6) 



Here we omitted the term with a full derivative propor- 
tional to NS(j> because it does not create the interference 
effects for the tunneling between opposite points of the 
unit sphere for a fully compensated antiferromagnet. 

At zero magnetic field instanton solutions that give 
minima of the action (J2J with the boundary conditions 
Z(±oc) — > ±e 2 can be easily found. The above-mentioned 
conditions take the form 9(— 00) = 0, 9(+oo) = it with 
arbitrary values for <^(±oo). The instanton solutions cor- 
respond to the rotation of I in the symmetry planes of 
the system, i.e., <f> = nk/2, where k is integer, and 



cos 1 



= tanh(cjT), 



sm( 



cosh(o;r) ' 



(7) 



where a = ±1 is the polarization of the instanton. 
The "frequency" ZU is 2S\/JzB u /h for the instanton 
paths going through the medium axis cf> = 0, it and 
ZJ = 2S \J J z(B u + B p )/h for the paths going through the 
hard axis <j> = ±7r/2. The Euclidean action calculated on 
these instantons is simply Ae — hN(huj/Jz). Following 
the saddle-point approximation that corresponds to the 
dilute instanton gas we take into account only paths with 
a minimal real part of the Euclidean action, i.e. with the 
smallest 57 = 25 V JzB u /h. Thus, at zero magnetic field 
the only instanton pair with a = ±1 and the rotation 
of I in the plane xz contributes to the tunneling. The 
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presence of the magnetic field drastically changes the fea- 
tures of the spin tunneling between two different classical 
ground states in the antiferromagnetic particle. 

First of all we consider the influence of the magnetic 
field on the static classical properties of the antiferromag- 
netic particle. We will use the expression (@} and restrict 
ourselves to the cases when the field is parallel to one of 
the symmetry axes x, y, or z. In these cases the field sim- 
ply renormalizes the value of the appropriate anisotropy 
constants. 

The influence of the field on the ground state is essen- 
tial when the field is directed along the easy axis z. In this 
case the renormalized constant of the uniaxial anisotropy 
B U (H) can change its sign, B U {H) = B u (l - H 2 /H 2 ), 
where 

H u = (8) 

9VB 

is the field of the familiar spin-flop phase transition in the 
classical theory of antiferromagnetism. When H > H u 
the states with I \\ ±e z become unstable and the ground 
states have I parallel to ±e x . 

When H is parallel to the medium axis x, the constant 
B u increases, B U (H) — B u + (g^ B H) 2 /(AJzS 2 ) and B p 
decreases as B p (H) = B p — (g^i B H) 2 /(AJzS 2 ), with the 
growth of the field. Thus, the field does not effect on 
the ground state, but a strong enough field can change 
the type of the axes x and y in the basal plane. When 
H > Hp, where 

g^B 

the axis y becomes an easy direction in the basal plane 
(medium axis) and the axis a; is a hard one. 

Finally, if the magnetic field is directed along the hard 
axis y, both anisotropy constants increase as B UiP (H) = 
B u .p + (g^sH) 2 / (4JzS 2 ). In this case, the field does not 
change the ground state, as well as the type of the axes. 

The naive substitution of the renormalized constants 
into the expressions Q) for an instanton and the Eu- 



clidean action calculated on it leads to a wrong predic- 
tion that when the field is directed along the easy axis 
Ae — ► at H — > H u , and for the field along the medium 
axis the values of Ae on the two classes of trajectories 
[<p = ixk and <\> = (2k + 1)tt/2) become equal at H = H p . 
As we will show below, both suggestions are wrong. 

In addition to the abovementioned renormalizations of 
the anisotropy constants, the field changes the dynamics 
of the vector I and leads to the appearance of a gyro- 
scopic term linear in dl/dr in the Lagrangian. The role 
of gyroscopic terms for tunneling in antifcrromagnets has 
been discussed in Refsi 12 ' 1 ? . The authors of these papers 
have shown that the gyroscopic term caused by the mag- 
netic field can create the imaginary part of the Euclidean 
action Ae that is a linear function of the magnetic field. 
As we have shown, similar effects can also arise due to the 
Dzyaloshinskii-Moriya interaction^ The imaginary part 
of the Euclidean action Ae can lead to the interference 
effects and to the oscillations of the tunnel splitting as a 
function of the magnetic field. We will show below that 
the gyroscopic term produces a dynamic renormalization 
of the real part of Ae, which is quadratic in the magnetic 
field. It can completely suppress the static contribution 
to Ae coming from the renormalization of the anisotropy 
constants. 



II. INSTANTON SOLUTIONS 

In order to describe macroscopic quantum tunneling 
between two classical states I — e z and I = —e z in the 
saddle-point approximation, it is necessary to find instan- 
ton solutions of the two Euler-Lagrange equations for the 
Euclidean action (J2J) for the independent variables 9(t) 
and 4>(t). Using the identities (5/59) J drH ■ (I x I) = 
-2<j>(Hl)sm6 and (6 /5(f)) J drH ■ (I x /) = 26 (HI) sm9 
that can be obtained through the variation of the La- 
grangian JfjJ), it is convenient to write down the Euler- 
Lagrange equations with the magnetic field H directed 
along an arbitrary symmetry axis 



- {lu 2 u + <j) 2 + j 2 (H 2 - H 2 ) + [uj 2 p + ~f 2 (H 2 - Hi)] sin 2 </»} sin(9cos(9 = -2ijr/{Hl) wa.6 , (10a) 
</> sin 2 9 + 24>9 sin 6 cos 6 - [uo 2 p + j 2 (H 2 - H 2 )] sin 2 9 sin <j> cos = 2i9~/(Hl) sin 6 . (10b) 

I 



Here u> u = jH u = 2S\ / 2J zB u /h and oj p — jH p = 
2Syj2 JzB p /h. The terms in the right-hand side are re- 
sponsible for the gyroscopic dynamics of the vector I. 
Note that we consider the field directed along one of the 
crystalline axes only, so only one of the components H 
in the system l|10[l is nonzero. 



It is important to note that in the case H ^ simple 
planar solutions such as <j> = irk/2 may not exist. If such 
solutions are absent, the full system (|10(l is equivalent 
to the Lagrange equations for a mechanical system with 
two degrees of freedom. To integrate such a system, the 
existence of two independent integrals of motion is nec- 
essary. For the general case lu u ^ and uj p ^ only one 
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first integral is known 
£ 



e 2 



b 2 sin 2 I 



u 2 sin 2 I 



2 • 2 n ■ 2 

UJp sin 9 sin 



0, (11) 



with the value £ = for separatrix solutions we are in- 
terested in. For this reason the system HlOf) cannot be 
solved analytically. However, approximate solutions can 
be constructed for all cases of interest; see a detailed con- 
sideration in the following sections. 



A. Field parallel to the easy axis 

We will start from the case of the field parallel to the 
easy axis, for which the gyroscopic terms in the right- 
hand side of Eqs. l(TU|l are independent of (j). The field 
does not violate the rotational symmetry around the easy 
axis, and the model with an isotropic basal plane (lo p — 
0) has a physical meaning for this case. Its analysis leads 
to instructive results and we give it completely. 

If uip = 0, the system (|10|) has one more integral of 
motion, which can be written as 



il = (cb — ijH) sin 9 = const 



(12) 



Using the two integrals of motion (|llf> and (|12fl . we can 
simplify the system to the ordinary differential equation 
for 8(t) only 



2 • 2 

to,, sin i 



O 2 



sin 2 9 



(13) 



and integrate it exactly. The instanton solutions with 
appropriate boundary conditions 9 — > 0, it and 9^0 
at r — * ±oo correspond to the values of the integrals 
Q, = 0, £ = 0. Thus, the solution is <f> — ijH(r — T\), 
where t\ is an arbitrary constant. Eq. (|13fl simplifies 
to 9 — uj 2 sin 2 9, and its solution is described by the 
formula with 57 = uj u . In terms of the vector I, the 
instanton solution for H parallel to the easy axis and 



takes the form 



l nr. 



cosh[7i7(r — ti)] 
cosh[ZZT (T - r )] 
sinh[7_ff (t — Ti)] 
cosh[cj (T - T )J 
l z = tanh[a7 (r - r )] . 

This solution has correct instanton asymptotics I 



(14a) 

(14b) 
(14c) 





and l z — > ±1 at r — > ±oo for all H in the range < lu u , 
which corresponds to H < H u . It is valid in the full 
region of stability of the states I — ±e z . Though the 
constant of effective anisotropy B U (H) changes from B u 
at H = to at H = H u , the Euclidean action calculated 
on the solution (|14|) does not depend on H, its real part 
is the same as at H = and the imaginary part is zero. 
Thus, in the uniaxial case B p = and H < H u the action 
is 



A E = 2hNS\ 



(15) 



This simple example shows that the gyroscopic term 
can essentially suppress the static renormalization of the 
anisotropy constant. For the case oj p = the renormal- 
ization is completely compensated by the gyroscopic term 
and Ae is real and does not depend on the field. 

For the case ui p = an exact solution of the problem of 
small fluctuations around the classical instanton solution 
is possible to construct, and the preexponential factor 
can be explicitly written. To do so, we introduce small 
perturbations z? and fi as 



(t)=0 o (t)+0, 
(r) = Mr) + 



sin O (t) 



(16a) 
(16b) 



where 9q(t) and </>o( r ) correspond to the instanton so- 
lution given by Eqs. (|14f> . The variational part of the 



action is Ae — Ag , where A K ^' is calculated on the un- 
perturbed instanton solution. It can be written as a sum 



(0) 



(0) 



5 2 A 



E 



of two independent terms Ae — A 
with the decoupled degrees of freedom $ and /j,. 
terms have the same form 



S 2 A 



(0) 



5 2 A 



{a) _ Nh z uJ 
AJz 



d{(,f {a) ,Mf {a) ) 



(ah 



Both 



(17) 



where £ = lot, is ^ or /i for 5 2 A% ] or S 2 A^\ 
tively, and M is the linear operator 

"de + "c^ii 2 ^' 



M 



respec- 



(18) 



This operator frequently appears in scattering problems 
associated with soliton theory, so its properties are well 
studied. The full set of its eigenvalues and normalized 
eigenfunctions is 



Mfo = 0, 



Mf k = (l + k 2 )f k , f k 



y/2 cosh £ 
(tanh£- ik)e~ tk ^ 



(19a) 



(19b) 



The mode /o is a localized eigenfunction and the modes 
with f = fk form a continuous spectrum. 

It is worthwhile to note that in contrast to other prob- 
lems of macroscopic quantum tunneling, for the case of 
antifcrromagnets with uj p = the zero mode /o is present 
for both kinds of fluctuations. The second zero mode 
is caused by the exact rotational symmetry around the 
axis z. The full preexponential factor for this prob- 
lem is a square of the usual fluctuation determinant 

D = y A^ /{2nK), and the tunnel splitting of the lowest 
levels A takes the form 



A = Chuj v 



A 



2nh 



exp 



(-Af/n) 



(20) 



where C is a numerical constant of order of unity. The ad- 
ditional large factor J A^ /(2nK) is a consequence that 
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the instanton solution (|14f> in the case uj p = contains 
two (not one, as usually) continuous parameters r\ and 
To. Note that the level splitting A does not depend on 
the magnetic field even if the fluctuation determinant 
is taken into account. This result can be explained us- 
ing exact quantum-mechanical arguments. For uniaxial 
system the Hamiltonian of the system commutes with 
the ^projection of the total spin S z , and the eigenstates 
of the problem can be characterized by definite values 
5 a = 0,±l,±2,... The two lowest levels form a doublet 
with S z = and the magnetic field does not influence 
on them. For this model the spin-flop transition at the 
field H — H u corresponds to the change of the value of 
S z for the lowest level from S z — to the value S z — 1 
or higher. Then, the growth of the magnetic field leads 
to the growth of S z and a sawlike dependence A (if) ap- 
pears. Since these effects are not associated with tun- 
neling, we will not consider them in the following, and 
restrict ourselves only to the region H < H u . 

Now we consider a more general case io v ^ 0. It 
is clear that at finite nonzero uj v the additional factor 
(Ae /27rft) 1 / 2 is absent, and the preexponential factor is 
smaller than that for the uniaxial case uj p = 0. The 
levels cannot be characterized by the quantum number 
S z = 0, ±1, ±2, and their splitting becomes depen- 
dent on the magnetic field. 

The instanton solutions with 9 = 0(t), <fi = exist at 
H = 0. It can be expected that at jH <c uj u the value 
of <j) or, more accurately, the appropriate projection of 
the vector I, l y ~ sin</>sin#o is small, and the function 
0(t) can be given by the solution similar to Eq. I|14|) . If 
ujp 3> uj u , the out-of-plane components of I are small, 
i.e. \l y \ <C 1. Assuming that \l y \ <C 1 for all values 
of parameters, the variational technique can be applied 
for evaluating the action. Choosing as a trial function 
cos 9 = tanh[rj(r — ti)] with some parameter uj that have 
to be found later, we rewrite Eq. IjlObj! in the linear ap- 
proximation in l y as 



trial parameter uj takes the form 



(M + e) l y (0 



2i~/H sinh £ 



cosh 2 £ 



(21) 



Here e = (uj p /uj) 2 , £ = ujt, and M is the linear opera- 
tor l|18|) . introduced above, with non- negative eigenval- 
ues. Since e > 0, the operator M + e is positively defined. 
It has an inverse operator which can be written through 
the ordinary bra and ket notation of eigenfunctions as 



1 



1/oX/ol 



M - 



St 



\fk)(fk\ 
+ e + k 2 



(22) 



where the summation is extended over the continuous 
spectrum. Since (sinh£/ cosh 2 £, / ) = 0, the formal so- 
lution of Eq. H21|) does not contain a term with 1/e. The 
solution Zj,(£) is determined by the summation over the 
states of the continuous spectrum only. After sim- 
ple calculations the Euclidean action as a function of the 



h 2 N 
TJz 

2uJ 3 



57 



dk 



(1 + fc 2 + uj 2 /uj 2 ) cosh^(7rfc/2) 



(23) 



If ojp — 0, the minimum of the action Il2i'>[l is reached 
at ZJ = uj u and we return to Eq. (|15|) again. Thus, in 
agreement with the previous analysis the Euclidean ac- 
tion depends on the field for uj p ^ only. If the ratio 
uj.pl uj u is small, this dependence is weak, and for the case 
of the field directed along the easy axis the Euclidean 
action is 



Ae(EA) = 2hNS\ 



<lH) 2 ujp- 



2uji 



(24) 



Thus, for ojp <§; oj u the dependence of the Euclidean ac- 
tion Ae on the field is weaker than it can be obtained 
from the naive consideration (see Fig. . The problem 
also can be solved in the limit case uj p 3> uj u , when the 
approximate solution of Eq. 121fl can be written as 



1 1, 



2i"fHoj sinh£ 
w 2 cosh 2 £ 



< 1 



(25) 



The field dependence of the Eu clidean action can be eval- 
uated as Ae - B U (H) = yj\ - {H/H u ) 2 . Thus, the 
naive consideration of the magnetic field through the 
renormalization of the anisotropy is recovered only in the 
limit case of high planar anisotropy uj p /uj u — > oo. 



B. Field perpendicular to the easy axis 

The case of the field directed along the hard axis is 
the simplest one. The plane xz remains preferable for 
the rotation of the vector I for all values of the field. 
The right-hand sides of the system l|l(Jfl are proportional 
to sin 2 9 sin 0, and the exact solution of the system 11U|) 
is <j> = 4>q = irk, uj = uj u . It corresponds to the rota- 
tion in the most preferable plane xz. The real part of 
the Euclidean action is independent of the field, but the 
imaginary parts have opposite signs for the two equiva- 
lent instanton trajectories with = and <j) — n, 

Ae(HA) = 2hNS^±i 7rH9fi f N . (26) 
V Jz 2Jz 

The nonzero imaginary part of the action H26[> leads to in- 
terference of the instanton trajectories and an oscillating 
dependence of the transition probability on the fieldii 2 ^ 
In this geometry the most interesting effect arises due to 
the fluctuation determinant. 13 ' 15 The equations for small 
fluctuations and fi around the instanton solution are 
uncoupled again, but a complex- valued potential appears 
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H c H„ Hsf H 



FIG. 1: The real part of the Euclidean action in units of 
h 2 N/2Jz as a function of the magnetic field (schematically). 
Data for the field directed along the easy, medium, and hard 
axes are marked as EA, MA and HA, respectively, close to 
the appropriate curves. Solid lines represent the calculations 
of Ae on the base of the full system HUH , dashed lines are 
results of the naive consideration. 



in the equation for The eigenvalues and the fluctua- 
tion determinant are also complex. Thus, the presence 
of Dg changes not only the tunneling amplitude for the 
instanton solution with a given value of cr, but also a 
phase shift for two instanton paths with a = ±1. At 
low fields Dg is almost real. It creates a decrease of the 
tunnel amplitude with the growth of the field. When the 
field increases, the factor Dg produces the phase shift of 
oscillations of A(H) caused by the interference. The au- 
thors of Ref^ suggest that such a complicated behavior 
of A(i?) caused by the fluctuation determinant Dg and 
the interference of instanton trajectories with ItiiAe 
can be a universal feature for a wide class of antiferro- 
magnets. As we saw earlier, for the field directed along 
the easy axis the dependence A(iJ) is strongly different 
from the scheme proposed by the authors of Refii^. For 
the field directed along the medium axis the behavior is 
also essentially different. 

The case of the field parallel to the medium axis is the 
most complicated one, but it is interesting. The right- 
hand sides of Eqs. (|10() are proportional to sin 2 #cos</>. 
One type of exact solutions with the rotation of I in 
the plane xy containing the hard axis y (a planar so- 
lution) can be written as <p = <po — Tt(k + 1/2) and 

cos 6 — tanh(cJr) with ZJ — i/wjj + ■ For this planar 

solution the structure of the Euclidean action Ag (MA) 
is the same as for the case of the field directed along 
the hard axis. Particularly, the real part of A^(MA) 
does not depend on the field and the imaginary part is 
proportional to the field 



A% ] (MA) = 2HNS 



B u + B, 

77 



p .irhgn B HN 
~ ±l 2Jz 



For such instanton trajectories the fluctuation determi- 
nant is similar to the factor obtained for the case of the 
field directed along the hard axis. It contains an imagi- 
nary part which leads to an extra contribution to the in- 
terference effects. But in contrast to the case of the field 
directed along the hard axis here the plane tfi — 7r(fc+l/2) 
is not the plane passing through the medium axis. The 

solution with <j) = (f>Q = ir(k + 1/2) and ZJ = \J + w p 

satisfies the system (|1()H exactly, but the real part of the 
action is not minimal, at least at small magnetic fields 
H < H p . In order to explain this fact it is sufficient to 
consider the value H = 0, when we have an exact so- 
lution with the rotation in the plane zx cf> — irk and a 
smaller real part of the action. Thus, only for high fields 
the planar solution can be relevant. For the naive consid- 
eration of the field the real part of the Euclidean action 
for the solution with <j) = ixk and 6 = 9(t) has to be pro- 
portional to ^ + (-yH) 2 , and Re A^ (MA) is equal to 

Re A^ (HA) at the point H = H p . But as we will see 
below, the situation is actually more complicated. 

Due to nonzero gyroscopic terms in the system l|10(l 
an exact solution with <j) — irk does not exist for the 
case H x ^ and an appropriate approximate instan- 
ton solution is nonplanar. Following the treatment of 
the preceeding section, we write the solution as cos 9 = 
tanh(cJT), ly = sin (j) sin 6 <^ 1 and determine l y from the 
linear equation 



(M + e) l y (0 



2ijH 



1 



cosh £ 



(28) 



where e = [cu 2 — (^H) 2 ]/uj 2 , £ = ujt, and M is defined 
by Eq. I|18fl . In contrast to the similar equation (|21|l the 
right-hand side of Eq. (|28|l is symmetric with respect to £ 
and a contribution from the localized eigenfunction l|19a|) 
is present. This contribution is proportional to 1/e and 
is mostly important for the case jH, uo p <C uj u . As we 
will see below, such an instanton is important at low field 
and the abovementioned restriction is irrelevant. The Eu- 
clidean action for this nonplanar instanton as a function 
of ZJ at ujp < uj u can be written as 



A E (to) 



K 2 N 



~f 2 H 2 



2lo 



(tt-jH) 2 to 



2(uf } 



7 2 iJ 2 ) 



AA E , (29) 



where AAe is a contribution from the continuous spec- 
trum that is determined through an integral over k with 
a structure that similar to Eq. (|23() . This term contains 
the factor ("jH/Zo) 2 and it can be omitted for the case 
of interest jH ~ uj p <g; Z3. After minimizing over ZJ the 
action for a nonplanar instanton for the case u> p < uj u 
takes the form 



A ( e P) (MA) = 2hNS 



(27) 
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where numerical constant r/ = ir 2 /4. The problem can be 
solved in the opposite limit case, namely, at u) u < u p . For 

this case one can replace the operator M + e by e, and the 
approximate solution reads l y (£) — 2ijH/(eLJcosh 2 £). 
Then, for the Euclidean action we arrive to the same 
equation (|30|) . but with another numerical constant r/ = 
8/3 w 2.667. Comparing this value found for lo u < u> p 
with 7r 2 /4 w 2.467 for w p < w„, we can tell that for the 

two opposite limit cases the Euclidean action A^ (MA) 



is approximately described by the same equation. Thus, 
we can suggest the Eq. (|3*5|) is a good approximation for 
any relation between u> p and ui u , that is in line with nu- 
merical data, see Sec. IIIII 

Thus, when the field is parallel to the medium axis, 
the Euclidean action is real for a nonplanar instanton, 
and due to Eq. (|30() it increases faster than it can be ex- 
pected from the static renormalization of the anisotropy 
constant. It is possible to show that both fluctuation 
determinants Dg and are also real. The interference 
effects are absent, and the tunnel splitting monotonically 
decreases with the growth of the field at small H < H Cl 
where H c is a critical field, at which the values of the Eu- 
clidean action for the nonplanar instanton becomes equal 
to the real part of A E 7 '(MA) for planar instantons. For 
small anisotropy in the basal plane, u p < w M , the value 
of H c is small, and 



H c — Hp 



B„ 



B p + (n/2) 2 B u 



< Hp 



(31) 



But it is smaller that H p even in the opposite limit case 

LO u < Up, 



H r — H„ 



1 



1 + 



0.6163R 



p 



(32) 



Thus, at low magnetic fields H < H c the Euclidean ac- 
tion for the nonplanar instanton solution A^ p ^ (MA) is 
lower than for the planar solution. Its value reaches 
A ( i ] (MA) at H = H c , and the scenario of tunneling is 
changed to the planar one, common to that is present for 
the field along the hard axis with the tunneling exponent 
independent of the magnetic field and with interference 
effects caused by imaginary parts of both Ae and fluc- 
tuation determinant. 



III. NUMERICAL DATA 

The semiclassical analysis of the coherent quantum 
tunneling between the classically degenerated ground 
states demonstrates that the level splitting is highly sen- 
sitive to the orientation of the magnetic field. Among 
the considered field orientations along the axes of rhom- 
bic symmetry the cases of the easy and medium axis are 
the most interesting ones. In both cases the Euclidean 
action has a zero imaginary part, and the corresponding 
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9^bH 
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FIG. 2: Splitting of the lowest level for the quantum model 
with Hamiltonian spin S = 5 and B u /J = 0.1 for three 
values of the ratio B v /B u = 0.0 (solid lines), 0.2 (circles), 
0.4 (triangles). The magnetic field is directed along the easy 
axis and normalized to the exchange field H ex = JzSN/(gfis) 
with 2 = 1 and N = 2. The cusps on the curves correspond to 
the change of the ground state that is a quantum counterpart 
of the spin-flop transition. 



interference effects are absent. We also have shown that 
for these two cases the prccxponcntial factor, which could 
be a source of interference, is real. Thus, for these orien- 
tations of the magnetic field interference effects does not 
appear, and the level splitting is mainly determined by 
the dependence of the real part of the Euclidean action 
on the magnetic field H. The character of this function 
is determined by the anisotropy B p in the basal plane. 
The dependence is absent for the case B p = and the 
field parallel to the easy axis only. The exponential fac- 
tor exp(— Ae/K) is an increasing function of H for the 
field parallel to easy axis, and a decreasing function of H 
for the field parallel to the medium axis. This behavior 
strongly differs from that is present for the field directed 
parallel to the hard axis. 

In order to check the semiclassical results found in 
assumption of some inequalities such as B p <C B u or 
B u <C Bp and to estimate the role of the fluctuation 
determinant, which was not investigated here, we diago- 
nalize numerically 30 the Hamiltonian Q for the two-spin 
quantum model with high enough values of the spin (up 
to S = 100), a small uniaxial anisotropy B u /J — 0.01 - 
0.1 that guarantees reasonable values of the level split- 
ting, and for several values of B p /B u , see figures below 
in this section. The numerical results show a satisfac- 
tory agreement with the instanton approximation even 
for moderate values of the spin S = 10 - 20 even with- 
out taking into account the preexponcntial factor. But 
some discrepancies, which cannot be attributed to the 
approximations used in the analytical consideration, are 
also seen. 

For the case of the field directed along the easy axis 
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FIG. 3: Splitting of the lowest level for the quantum model 
with Hamiltonian JQ, sphi S — 5, and B U /.J — 0.1 for three 
values of the ratio B p /B u — 0.0 (solid lines), 0.2 (circles), 
0.4 (triangles), and the magnetic field directed along the hard 
axis (open symbols) or medium axis (full symbols). The field 
is normalized to the exchange value. 



at Bp = the tunnel splitting is independent of the field 
up to the point of the spin-flop transition. This behavior 
is exactly reproduced by means of a numerical diagonal- 
ization, see Fig. [3 Due to the results of the instanton 
approach, for B p ^ the Euclidean action decreases, 
and the level splitting increases with the growth of the 
field H . The effect becomes more pronounced at large 
Bp. This behavior well corresponds to the numerical 
data represented in Fig. The only difference between 
the numerical data for the level splitting E\ — Eq and 
the semiclassical results is that the field of the spin-flop 
transition identified as a field of the cusp on curves in 
Fig. |3 slowly depends on the ratio B p /B u for the quan- 
tum model. This dependence is completely absent in the 
semiclassical approximation, in particular, in the instan- 
ton approach. Probably, close the point of the spin-flop 
transition, where the sharp decrease of the Euclidean ac- 
tion is present, the value of Ae/^ becomes comparable 
with unity even for large S, and the quantum fluctuations 
treated beyond the semiclassical approximation become 
important. 

In the case of the field directed along the hard axis 
as well as in the case of high fields, H > H c , along 
the medium axis the instanton approach predicts that 
the tunneling occurs through the planar instanton paths. 
For these cases the level splitting A(ff) oscillates as a 
function of the field with a constant period. If the field 
is directed along the hard axis, the behavior of A(i7) 
coincides with the results of Refsii^i^. The amplitude 
of oscillations is determined by the preexponential factor 




H/H p 

FIG. 4: Numerical data for the logarithm of the level split- 
ting normalized by the value at H = versus the magnetic 
field directed along the medium axis for the spin S — 10 and 
B u /J = 0.1 and three values of the anisotropy constant in 
the basal plane: B v /J — 0.05 (circles), 0.15 (squares), 0.25 
(triangles). The magnetic field is normalized to the reori- 
entation field H p . The solid line describes the theoretical 
dependence (1301) for a nonplanar instanton. 



only, and it weakly depends on B p . This kind of behav- 
ior is clearly seen in Fig. [3] In the case of high fields, 
H > H c , directed along the medium axis the tunneling 
is also determined by planar instantons with a nonzero 
imaginary part of the action and interference effects ap- 
pear. But when the ratio B p / B u increases, the main ex- 
ponential factor drastically decreases, and the amplitude 
of oscillations A(iJ) decreases also. This feature is in a 
good agreement with the numerical calculations depicted 
in Fig. 

At low fields, H < H c , parallel to the medium axis the 
tunneling is determined by nonplanar instantons. The 
real part of the Euclidean action Re A^ p ^ [MA) mono- 
tonically increases, and the instanton approach predicts 
a strong monotonic decrease of A(H) up to H = H c . For 
any values of B p /B u at H < H c the instanton approach 
also predicts that oscillations caused by interference are 
absent. The effect of the axes reorientation in the basal 
plane at H — H p does not clearly manifest itself as the 
spin- flop transition in Fig. [21 but it can be understood as 
a shift of a point, where oscillations start, to the region 
of high fields. Both factors, as well as the growth of the 
characteristic field of transition to the high-field tunnel- 
ing picture H c , are in a qualitative agreement with data 
of numerical calculations presented in Fig. [3] 

In order to give a more detail comparison of the an- 
alytical and numerical data for small fields, H < H Cl 
parallel to the medium axis, we investigated numerically 
the value of spin S = 10, for which the role of the field 
dependence of the fluctuation determinant is expected to 
be less important. The data together with the simple 
theoretical estimate of the level splitting, see Eq. (|3(J[1 , in 
the form A(if)/A(0) = exp{-[A(H) - A(0)]/H}, where 
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FIG. 5: Real part of the Euclidean action per unit spin 
obtained from the quantum model with Hamiltonian Q, 
B u /J — 0.1 for three values of the ratio B p / J = 0.02 (cir- 
cles), 0.04 (squares), 0.08 (triangles) and the magnetic field 
directed along the easy axis. Solid lines are analytical pre- 
dictions, Eq. 12411 . The field is normalized to the exchange 
value. 
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FIG. 6: Real part of the Euclidean action per unit spin 
obtained from the quantum model with Hamiltonian (QQ, 
B u /J = 0.1 for three values of the ratio B v / J = 0.05 (circles), 
0.1 (squares), 0.2 (triangles) and the magnetic field directed 
along the medium axis. The magnetic field is normalized to 
H p for each value of B p . Solid lines are analytical predictions, 
Eq. 13011 with H p as a trial parameter. 



A(H) is the value of the Euclidean action for nonplanar 
instantons, are present in Fig. Here we can say about 
at least a semiquantitative agreement. 

In order to check quantitatively the analytical expres- 
sions (|24[l and l|30l) for the real part of the Euclidean 
action, we propose the method of extracting the value of 
Ae from the level splitting data without using an explicit 
expression for the preexponential factor. This method 
is based on the general theoretical formula for the level 



splitting 

A = DS f3 exp(-SA) , (33) 

where D, (3 and A are functions independent of the spin 
with normalized values of the magnetic field H/H p or 
H/H u . In the semiclassical approximation the function 
A means the real part of the action per unit of spin, D 
is the preexponential factor and (3 is associated with the 
number of zero modes: (3 = 2 for the case of the purely 
uniaxial symmetry, that is B p = and the magnetic field 
parallel to the easy axis, see also Eq. l|20[l and discussion 
there; (3 = 3/2 for all rest models with biaxial symmetry. 

The parameter (3 is also chosen floating in order to 
absorb corrections to the cr-model and its semiclassical 
treatment. This method could, in principle, be applied 
for any semiclassical problem, but it is mostly useful for 
problems, where the imaginary part of the action, as well 
as the preexponential factor, are zero, and interference 
effects with oscillations are absent. For the tunneling in 
the external field this is just the problem we are interested 
in. This condition is realized at H < H u for the field 
directed along the easy axis and H < H p for the case of 
the field parallel to the medium one. 

In order to obtain the values of interest, a numerical 
calculation of splitting of the lowest level is performed for 
few (at least, three) large values of spin. The real part 
of the action can be easily extracted from the numerical 
data such as qualitatively presented in Figs. [S] and The 
functions Ae{H) in these figures are obtained by fitting 
the numerical data for the spin 10 < S < 20 and each 
value of the field using Eq. . 

For the simplest uniaxial case B p = and the field 
directed along the easy axis we found that the expected 
value (3 = 2 appears, and the Euclidean action is in- 
dependent of the field up to the point of the spin-flop 
transition. Then, with growth of B p we found a very 
sharp transition to the regime with one zero mode, for 
which the value of (3 = 3/2 is reproduced. Thus, checking 
that (3 does not differs significantly from 3/2 (the range 
1.5 ± 0.1 is taken) side by side with more trivial condi- 
tions that the obtained values of D, (3 and A are stable 
with respect to varying S and the exponent SA is large 
enough, we select sets of the parameters for the quantum 
model (UJ). 

The action A obtained by fitting of Eq. l.'i.'ill with the 
fixed value (3 = 1.5 is plotted as functions A(H) in Figs. 03 
and El for the fields directed along the easy and medium 
axes. Note that the action A(0) at zero field is sub- 
tracted from the functions A(H). The analytical theory 
gives the value A(0) = A^/B u / J independent on B p , but 
numerical calculations demonstrate a weak dependence 
of A(0) on B p . In details, A th (0) = 1.26 for B u /J = 0.1 
and A num (Q) is 1.22, 1.19 and 1.14 for the same value 
B u /J = 0.1 and B p / J equal to 0.05, 0.10, 0.20, respec- 
tively. 

The only way to describe the numerical data is to con- 
sider the fields H u and H p in Eqs. I|24|l and (|30|l as phe- 
nomenological parameters that predicted by the classical 
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expressions (jSJ) and For simplicity, we rescale the 

fields in Fig. as H — > H/H p . The numerical data are 
fitted by Eqs. 1)24(1 and 1)30(1 using these fields in the range 
H < 0.1H U and H < 0.2H p . Obtained values for a trial 
parameter H p /H p = 1.164, 1.144, 1.145 for B p /J = 0.05, 
0.10, 0.20 are in a good agreement with the classical re- 
sult, where obviously H p = H p . Appropriate analytical 
curves are plotted as solid lines in Figs. El and El In the 
case of the medium axis the theory predicts the shape 
of the curves obtained from the quantum model with a 
good accuracy up to fields close to H c . For the easy axis 
the action decreases more significantly that it would be 
expected from the perturbative treatment of the semi- 
classical model. In both cases we can pretend on the 
quantitative agreement of the proposed theory and nu- 
merical data. 

It is important to note one more discrepancy between 
the developed analytical theory and the presented numer- 
ical data. The analytical expression for the level splitting 
in the field directed along the hard axis does not contain 
any dependency on B p , but in Fig. 01 as well as in the nu- 
merical data of Refsi 13 i 15 , this dependency is present. It 
is more important for higher values of B p /B u . To explain 
it, as well as the observed dependence of the parameters 
H p and H u , we note that the tr-model treats antiferro- 
magnets in the first approximation over small ratios of 
the anisotropy constants or the magnetic field to the ex- 
change integral J. Here, the values of these ratios was 
taken in the range 0.1 - 0.2, and the deviation of the o- 
model results that is of order of 10% from the numerical 
calculations are not surprising. 

IV. CONCLUDING REMARKS 

In conclusion, the antiferromagnetic particles can show 
a reach variety of tunneling behaviors that depend on 
the direction of the magnetic field. In addition to the 
oscillation behavior for the field directed along the hard 
axis^^ we found the growth of the tunnel splitting A 
for the field directed parallel to the easy axis and a steep 
decrease of A for the field along the medium axis. Both 
mentioned behaviors are connected to the tunnel expo- 
nent dependence on the field. It is important to note that 
such effects cannot be directly associated with the de- 
crease or increase of the tunnel barrier, respectively, that 
governed by the static renormalization of the anisotropy 
energy. 

Let us briefly discuss the possibility for experimental 
investigations of the tunneling effects predicted in the 
paper. The main point consists in kinds of antiferromag- 
nets that could be used for experiments. The traditional 
antiferromagnetic samples such as small ferritin particles 
have unpaired spins and behave as to noncompensated 
antifcrromagnets. For this reason the destructive inter- 
ference for them is mainly dictated by the excess spin 
in the way common to ferromagnets^ii Moreover, the 
presence of a nonzero total magnetic moment drastically 



changes the structure of the ground state. It is enough to 
say that the degeneracy is absent except some fixed direc- 
tions of the field with respect to the crystalline axis4Sii£ 
We proposed a way to overcome this problem^ but lim- 
itations caused by noncompensated spins seams to be 
more serious. Note that the same problem appears for 
ferromagnetic particles where the effects of the barrier 
reduction^ and the oscillation behavior of the ground- 
state tunnel splittingSii 2 ^ was predicted many years ago, 
but observed only recently^ 

The key point in this important experimental success 
is based on the synthesis of high-spin molecules packed in 
the well-oriented monocrystals. Up to our understand- 
ing, the first possibility to investigate purely antiferro- 
magnetic features is to use high-spin molecules with a 
well-defined spin structure. The molecules with ferro- 
magnetic and antiferromagnetic couplings, uniaxial and 
rhombic anisotropics have been synthesized in the re- 
cent years^* For known ferromagnetic molecules such 
as Fes the splitting is small, but the technique devel- 
oped by Wernsdorfer and Sessolli^ allows one to mea- 
sure a very small tunnel splitting of order of 10 _8 K. The 
first possibility discussed by many authors consists in 
using spin rings with antiferromagnetic coupling. For 
well-known antiferromagnetic molecular magnets such 
as Feio, Fe6, Vs the problem is opposite to that for 
ferromagnetic molecules: the anisotropy is too small, 
and the barrier is too low to see clear semiclassical ef- 
fects such as MQTi 2 ^i 2 9 On the other hand, antifer- 
romagnetic rings of eight chromium ions with a high 
anisotropy have been recently synthesized^ One more 
possibility is to use spin dimers containing two coupled 
high-spin molecules (molecular magnets) with a ferro- 
magnetic coupling inside the molecule and an antifer- 
romagnetic intermolecular coupling. For instance, the 
observation of the well-structured dimers of high-spin 
molecules Mn4 (spin S = 9/2) with the antiferromagnetic 
coupling between two Mn4 molecules has been recently 
reported. 28 Quantum tunneling in the Mn4 dimers was 
investigated experimentally. 29 Such dimers of high-spin 
molecules such as Fes with the macroscopic spin S = 10 
and well pronounced rhombic anisotropy could be a good 
candidates for observation of the effects considered in our 
paper.' 51 It is worth to note also that the predicted pos- 
sibility of enlarging the value of the Euclidean action (to 
suppress the tunneling) by means of the magnetic field 
directed parallel to the medium axis can be useful for 
following investigations. 
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